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Abstract 

It is known that scintillation efficiency of organic single crystals depends on the 
direction of nuclear recoils relative to crystallographic axes. This property could be 
applied to the directional WIMP dark matter detector. 

The scintillation efficiency of carbon recoils in a stilbene crystal was measured 
for recoil energies of 30 keV to 1 MeV using neutrons from '^Li(p,n)'^Be and ^^^Cf. 
Anisotropic response was confirmed in low energy regions. The variation of the 
scintillation efficiency was about 7 %, that could detect the possible dark matter 
signal. 
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1 Introduction 



It is considered that the galactic halo is composed of weakly interacting mas- 
sive particles (WIMPs) as dark matter. These particles could be directly de- 
tected by measuring the nuclear recoils produced by their elastic scattering 
off nuclei in detectors [1]. However, nuclear recoils produced by background 
neutrons are indistinguishable from those by WIMPs. Therefore, one should 
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look for statistical signature of the WIMPs. Realistic distinctive WIMP sig- 
nals arise from the earth's motion in the galactic halo. Annual modulation 
of event rate is one of the possible WIMP signals caused by earth's revolu- 
tion around the Sun (~30 km/s). Actually, this signature was argued by the 
DAMA experiment [2] . 

Nevertheless, the most convincing signature of the WIMPs appears in the 
directions of nuclear recoils. It is provided by the earth's velocity through the 
galactic halo (~230 km/s). Assuming the WIMP halo is an isothermal sphere, 
the strong WIMP wind is blowing on the earth and the distribution of the 
nuclear recoil direction shows a large asymmetry. Hence, detectors sensitive 
to the direction of the recoil nucleus would have a great potential to identify 
WIMPs, and further, to provide information on the galactic halo [3]. 

It is known that scintillation efficiency of organic crystals to heavy charged 
particles depends on the direction of the particles with respect to the crys- 
tallographic axes [4,5]. This property makes it possible to propose a WIMP 
detector sensitive to the recoil direction of the nucleus [6,7]. The directional 
scintillation efficiency to recoil protons and carbons produced by neutrons in 
MeV regions was reported [8] , but the recoil energy given by WIMPs is much 
lower. 

We had measured the proton recoils in a stilbene crystal for recoil energies 
of 300 keV to 3 MeV and shown that the efficiency depends on the direction 
of the recoil proton [9] . However, carbon recoils would be more effective than 
proton recoils in detecting WIMPs with the interesting mass region of above 
10 GeV. In this paper, we report on the anisotropic scintillation efficiency 
of carbon recoils in a stilbene crystal with neutrons especially in low energy 
regions. We also estimate the sensitivity to the WIMP wind of the stilbene 
crystals. 



2 Experimental Setup 

The crystal lattice of stilbenes is shown in Fig. 1. Stilbene crystals form mon- 
oclinic systems and the crystallographic axes are called a, b, and c. The axis 
perpendicular to a-b plane is called c' . The direction of c' axis is easily dis- 
cerned because stilbene crystals are cleaved along a-b plane, and according to 
[5], the scintillation efficiency of stilbene crystals depends on the recoil angle 
with respect to c'. Therefore, we measured the recoil angle dependence by 
changing the angle of the recoil direction with respect to c' axis (9) and the 
angle around the c' axis {(p, (p — 0° was determined arbitrarily). 

In order to obtain high statistics and various incident neutron energy, two 
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neutron sources, ^Li(p,n)^Be and ^^^Cf were employed. The ^Li(p,n)^Bc source 
run was performed at 3.2 MV Pelletron accelerator of the Research Laboratory 
for Nuclear Reactors at Tokyo Institute of Technology. Pulsed proton beam 
interacted with a thin lithium target, and pulsed neutrons were produced by 
the ^Li(p,n)^Be reaction. The repetition rate of the pulsed proton beam was 
set at 2 MHz and the pulse width was 1.5ns. Changing proton energy three 
times, neutrons with energies of 200 to 650 keV were obtained. The beam and 
detector geometry is illustrated in Fig. 2. The temperature of the room was 
kept at 26°C during the measurement. 



The dimension of the target stilbene that we measured was 2 cm x 2 cm 
X 2 cm. Two opposite faces were cleaved {a-b plane) and other four sides 
(arbitrary plane) were polished. The polished sides were covered with GORE- 
TEX®, and 51 mm PMT (Hamamatsu H6411) was attached to each cleavage 
plane. Self- coincidence between two PMTs was required to reduce dark current 
events. Scattered neutron was detected by 51 mm 0x51 mm liquid scintillator 
encapsulated in an aluminum cell (Saint-Gobain BC501A-MAB1) with PMT 
(Hamamatsu H6411). The scattering angle was fixed at 120° to avoid proton 
recoil events. 



The incident and scattered neutron energies were measured by time-of-flight 
(TOE) method and recorded by a TDC (Hohshin C021). The coincident stil- 
bene output was used to define the common "START". The outputs of the 
BC501A detector and the delayed pulsed proton signal from the "Time Pick-off 
Unit" of the accelerator were used as the "STOP" signals. The PMT outputs 
were recorded by charge ADCs (Hohshin C009H) with two different gates — 
20 ns and 500 ns from the rise of the waveforms of the PMT outputs. Both 
stilbene and BC501A provide pulse shape discrimination (PSD) capabilities 
for n-7 separation arising from the difference of the slow scintillation compo- 
nent fraction [4] . Consequently, the ratio of charge integrated by two different 
gates gives n-7 information. 



The setup of the 252Qf source run was almost the same as the ^Li(p,n)'^Be 
source run, except the measurements of incident neutron energy. The tar- 
get stilbene crystal was placed 60 cm away from the 3.5 MBq ^^^Cf and the 
BC501A detector was placed 60 cm away from the target at 120° of the scat- 
tering angle. To determine the timing of the nuclear fission of ^^^Cf, prompt 
gamma rays produced by the fission were detected with the 2 cm thick plastic 
scintillator located by the neutron source. PMT (Hamamatsu R1250) outputs 
were delayed and used as the "STOP" signal of the TDC. 
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3 Measurement results 



The detected electron equivalent energy (visible energy) for carbon recoils in 
the stilbene crystal was calibrated with 5.9 keV, 6.4 keV, 14.4 keV, 22.2 keV, 
and 32.2 keV X/gamma rays from ^^Fc,^''Co/'^*^Cd, and ^^'''Cs. The efficiency 
close to the threshold was checked by the 5.9 keV X rays. We see from Fig. 3 
that the efficiency was not lost above 3keV. 

Although about 40% of neutrons interacts with the stilbene crystal more than 
twice in the measured energy region ^ , the single recoil events of neutrons 
scattered at 120° by carbons were selected using the TOFs of the incident and 
scattered neutrons. 

However, in the ^Li(p,n)^Bc source run, there still remained background events. 
These were gamma ray events from ^Li(p,7)^Be (minority) and nuclear recoils 
in the stilbene crystal accidentally coincident with the events in the BC501A 
detector (majority). The gamma ray events were rejected by means of PSD of 
the both stilbene and BC501A detectors. The events around ns of the TOF 
spectra of both incident neutrons and scattered neutrons are also rejected as 
the gamma ray events. On the other hand, to estimate the backgrounds of 
the accidental nuclear recoils, the measurements without requiring the coinci- 
dence of the BC501A detector were performed before and after each of the 6'/0 
run. Those background spectra were stable throughout the experiment. The 
stability of the beam intensity was also confirmed by the Li glass scintillator 
in the beam line. 

The visible energy spectrum with 9 — 90° for recoil energies of 100 to 105 
keV is shown in Fig. 4 together with the background spectrum after gamma 
ray events rejection. As indicated, background spectrum was normalized and 
subtracted. The normalization factor was the ratio of the number of events of 
the higher energy regions. The shape of the background spectrum below 10 
keV reflects the fact that the scintillation efficiency of proton recoils increases 
at low energies [10] (Not the loss in efficiency near the threshold). The visible 
energy of that recoil energies was derived by fitting Gaussian to the spectrum 
above 3keV. 

The analysis procedure of the ^^^Cf source run was detailed in [9]. After the 
events were binned with recoil energy, the visible energy was determined by 
taking the average in each bin of the recoil energy because of the small statis- 
tics. 

The derived scintillation efficiency relative to the efficiency for electrons with 
9 — 90° and ^ = 0° are shown in Fig. 5. The result is consistent with the mea- 

^ Based on the GEANT3 simulation 
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surements for proton recoils [9] and for charged particles in high energy regions 
[5], that is, the scintillation efficiency of stilbene crystal for carbon recoil is 
maximal in the direction perpendicular to the c' axis {6 = 90°) and minimal in 
the direction parallel to the c' axis {9 = 0°). On the other hand, ^-dependence 
of the scintillation efficiency was not observed in this measurement. That is 
also consistent with [5]. 

In order to estimate the response to the WIMPs, it is necessary to express 
the scintillation efficiency as a function of the recoil energy. The relation be- 
tween the light yield (dL/dx) and the energy loss (dE/dx) of scintilltors had 
been explained by Birks's empirical model [4]. However, this model does not 
describe the rise of the efficiency at low energies of our results. Recently, such 
enhancement at low energies for organic liquid scintillator was reported [10], 
and the light yield is introduced as 



with kBe/kBn = 3250. This relation is based on the assumption that electronic 
energy loss and nuclear energy loss contribute differently to the quenching 
process and scintillating process. The (dE/dx) terms are given by SRIM2003 
package [11] as a function of recoil energy. 

While Sg/ Sn, kBe, or both of them could be function of the angle 9 for this 
experiment, it is suggested that the response anisotropy is due to the difference 
in the value of the ionization quenching parameter kB [4] . Therefore, we fitted 
equation (1) to both of the measured 9 = 90° and 9 = 0° data at a time with 
the three parameters; kBe{9 = 90°), kBe{9 = 0°), and Se/Sn- The lines in Fig. 
5 are the best fits to the data with kBg{9 = 90°) = 14.3 ± 0.4 mg/cmVMeV, 
kBe{9 = 0°) = 15.3±0.5 mg/cmVMeV and SJSn = 0.547±0.030 (x^ = 68.7, 
d.o.f. = 43). Although statistics is not very high, anisotropy of 7% is seen over 
the measured energy region. 



4 Discussion 

The WIMP signal by the stilbene scintillator can be obtained by comparing 

the visible energy spectra measured for different orientation with respect to 
earth's motion in the galactic halo. The variation amplitude (difference of 
the spectra) should be independent of any terrestrial backgrounds. Assuming 
spherical isothermal halo model, we estimate the expected variation of the 




(1) 
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visible energy spcctrTim. The differential angular rate of WIMPs in detectors 
with respect to laboratory recoil angle 7 and recoil energy Er with WIMP 
mass and target nucleus mass M„ is given by [12]: 



,^ , = exp 

dERd cos 7 27r 2 M3 M„ 



{ve cos 7 



0^ 







(2) 



where po is the local WIMP density, Vq is the velocity dispersion of the isother- 
mal halo, and f^; is the component of the earth's velocity parallel to the galac- 
tic rotation; w^^^ = (M^ + M^f Eji/2M^Mn is the minimum WIMP velocity 
that can produce En and a is the WIMP-nucleus cross section. 

We calculate by Monte Carlo method in two cases, c' axis maintained to the 
direction of the galactic rotation (the direction of the constellation Cygnus) 
and to the galactic center (the direction of the constellation Sagittarius). 9- 
dependence of {dL/dx) is assumed to be trigonometric as observed in high 
energy region [5] . The expected spectra and the variation amplitude is shown 
in Fig. 6. This implies that it is plausible to use stilbene crystals for WIMP 
search. The variation of the spectra is much the same as that of annual mod- 
ulation. However, less systematic uncertainty of the measurements can be 
achieved because the direction of the crystal axis to the WIMP wind can be 
controlled. 

Other organic crystals, such as anthracene and naphthalene show stronger 
dependence of scintillation efficiency on direction to crystallographic axis [13]. 
These crystals should be more sensitive to WIMP wind and the variation of 
the spectra should be larger. Especially, it was reported that the variation of 
the scintillation efficiency of naphthalene crystal is about 50 %. With these 
organic crystals, it would also be possible to obtain some information on the 
motion of the galactic halo. 



5 Conclusion 



We measured the scintillation efficiency of carbon recoils in a stilbene crystal. 
The efficiency was observed to increase at low energies and depend on the 
recoil direction with respect to the d axis. The variation of the directional 
anisotropy is about 7 %. Assuming this anisotropic response, we estimate the 
sensitivity to the WIMPs and found that stilbene scintillator could detect 
robuster WIMP signal arising from the earth's rotation around the galactic 
center as compared with the annual modulation. 
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c' 



c 




Fig. 1. Schematic drawing of the stilbene crystal lattice and the definition of the 
recoil angles. 6 is the recoil angle with respect to the c' axis and ^ is the angle 
around the d axis. 
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Fig. 2. Setup of the ^Li(p,n)^Be source run. Because the crystal shape is cubical 

and c' is perpendicular to cleavage planes where PMTs are attached, the geometrical 
configuration of the crystal relative to the neutron beam is same for 6 = 90° and 
6 = 0°. 
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Fig. 4. (a) The visible energy spectra with 9 = 90° for recoil energies of 100 to 
105 keV (solid) and the normalized background (dashed), (b) The visible energy 
spectrum after background subtraction and the fitted Gaussian. 
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Fig. 5. The measured scintillation efficiency relative to the efficiency for electrons 
with 6 = 0° and 6 = 90° recoils = 0°). The lines show best fits to the data using 
equation (1). Se/Sn is 0.547 it 0.030 for the both cases. The inset is the results of 
the ^Li(p,n)'^Be source run. Horizontal error bars represent selected recoil energy 
region for calculating the efficiency. 
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Fig. 6. (a) The expected energy spectra of the stilbene scintiUator when d axis 
maintained to the direction of galactic center (sohd) and galactic rotation (dashed). 
The parameters that we used in the calculation were po = 0.3 GeV/cm^, v^^ = 220 
km/s, ve = 232 km/s, WIMP-proton spin independent cross section o-^-p = 10~^ 
pb, and = 50 GeV. (b) Difference in count rates and the asymmetry (= differ- 
ence/average). The experimental resolution is not taken into account in the plots. 
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